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SUMMARY 


A  theoretical  analysis  of  the  performance  of  a  propeller- 
duct  combination  operating  in  ground  effect  was  extended  from 
the  work  conducted  under  Contract  DA  44-177-TC-524,  which  was 
reported  in  Reference  1.  A  propeller  was  designed  for  optimum 
twist  and  taper  and  to  provide  a  constant  velocity  across  the 
GEM  nozzle,  thus  producing  a  maximum  lift-to-power  ratio  of 
the  GEM  in  accordance  with  this  reference.  The  optimum 
propeller  design  was  of  the  same  propeller  solidity  as  a 
rectangular,  non  twisted  blade  propeller  tested  previously 
to  obtain  a  direct  comparison  of  performance.  The  method  of 
analysis  used  to  determine  the  performance,  loads,  and  moments 
produced  on  a  blade  of  arbitrary  planform  and  twist  operating 
in  a  duct  is  presented. 

A  test  program  was  conducted  to  determine  the  effect  of 
the  optimized  propeller  on  GEM  performance.  The  duct  nozzle 
area  was  also  varied  to  determine  the  effect  of  this  parameter 
on  propeller  optimization. 

Propeller  blade  stall,  along  with  large  inlet  losses  which 
occurred  during  certain  operating  test  conditions,  resulted  in 
performance  which  could  not  be  accounted  for  in  the  analysis. 


1 


CONCLUSIONS 


The  thrust  produced  and  power  required  by  a  propeller 
operating  in  a  GEM  can  be  predicted  within  15 -percent  pro¬ 
vided  that  the  operating  condition  does  not  represent  condi¬ 
tions  for  which  the  propeller  is  stalled  or  the  inlet  losses 
become  large.  The  conditions  for  which  the  theory  does  not 
apply  overlap  the  normal  operating  conditions  of  some  GEMS, 
and  therefore  the  development  of  the  analysis  should  be 
continued  to  include  these  conditions. 

The  optimized  propeller,  which  produced  a  uniform  nozzle 
velocity,  showed  no  improvement  in  GEM  performance  over  the 
non  twisted  propeller  which  produced  a  triangular  nozzle 
velocity  distribution.  The  straight-bladed  propeller  tends 
to  give  better  inlet  performance.  The  internal  efficiency  of 
the  straight-bladed  propeller  was  10 -percent  greater  than 
the  optimized  propeller,  for  heights  above  h/d  of  0.10.  The 
optimized  propeller  produced  higher  power  loadings  and  inter¬ 
nal  efficiencies  than  the  non- twisted  propeller  below  h/d 
of  0.10.  This  improvement  is  principally  due  to  the  greater 
resistance  to  stalling  of  the  optimized  propeller. 

Duct  efficiencies  on  the  order  of  75-percent  were 
achieved  for  the  straight  duct  for  jet- to- total-base  area 
ratios  from  0.1  to  0.8.  The  duct  efficiency  is  dependent 
upon  the  total  pressure  of  the  flow  at  the  GEM  nozzle. 
However,  the  duct  loss  factor,  which  depends  upon  dynamic 
pressure,  increases  with  decreasing  jet  area.  The  duct 
losses  are  attributed  to  a  large  inlet  loss. 

The  configuration  in  which  the  jet  area  was  equal  to  60- 
percent  of  the  total  base  area  produced  the  highest  power 
loadings  of  the  configurations  tested.  This  configuration 
also  showed  the  best  performance  for  h/d  greater  than  0.1 
based  on  the  GEM  power  factor. 


2 


INTRODUCTION 


This  program  was  a  continuation  of  the  effort  being 
directed  toward  the  improvement  of  the  internal  efficiency 
of  GEMs.  The  development  of  the  GEM  has  demonstrated  a 
need  for  simple,  lightweight  propellers  and  ducts  which 
have  adequate  performance.  However,  there  have  been  in¬ 
sufficient  test  data  or  analysis  to  evaluate  duct  and  propeller 
designs.  The  effort  presented  in  this  report  is  directed 
toward  providing  a  guide  for  designers  in  optimizing  GEM 
propellers  and  duct  systems. 

Analysis  was  performed  to  develop  methods  of  calculating 
the  performance  of  an  arbitrary  GEM  propeller  as  well  as 
optimized  propeller  design.  Experiments  were  conducted  for 
comparison  with  the  analytical  treatment.  The  configuration 
tested  consisted  of  a  axially  symmetric  straight  duct  in 
which  the  jet  area  was  varied  by  the  addition  of  centerbodies 
of  various  diameters.  Two  propellers  were  tested  in  the  duct: 
one  propeller  had  a  constant  chord  and  zero  twist,  herein 
called  the  straight- bladed  propeller;  the  other  propeller 
tested  represented  the  optimized  propeller  designed  with  twist 
and  variable  chord  to  produce  a  uniform  axial  velocity  distri¬ 
bution  with  a  minimum  of  losses. 

This  report  is  divided  into  four  sections  and  three 
appendices.  In  Section  1,  the  development  of  propeller  per¬ 
formance  equations  formulated  during  this  program  is  presented 
along  with  the  analysis  for  the  determination  of  the  perfor¬ 
mance  of  a  propeller  with  an  arbitrary  planform  and  twist. 
Considerations  for  the  design  of  an  optimized  propeller  are 
discussed  in  Section  2.  The  results  of  the  experimental 
program  are  presented  in  Section  3.  A  comparison  of  these 
results  with  theory  is  given  in  Section  4.  The  test  pro¬ 
gram  and  a  description  of  the  instrumentation  are  presented 
in  Appendix  I.  In  Appendix  II,  an  analysis  for  the  deter¬ 
mination  of  loads  and  moments  produced  on  a  twisted  and 
tapered  propeller  blade  operating  in  a  duct  is  developed. 

The  design  of  the  optimized  propeller  used  in  this  test 
program  is  given  in  Appendix  III. 
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1.  PROPELLER  DESIGN 


1.1  Development  of  Propeller  Performance  Equations 


The  following  development  represents  a  method  of 
predicting  propeller- in-duct  performance  for  a  propeller  of 
arbitrary  geometry  at  any  operating  condition.  The  process 
also  results  in  a  predicted  axial  velocity  distribution. 

The  development  is  a  continuation  of  the  work  performed  in 
Reference  1.  The  theory  has  also  been  modified  to  include: 

a.  An  inlet  energy  loss  factor 

b.  Airfoil  section  performance  data,  thereby  including 
the  effects  of  Mach  number  and  a  small  amount  of 
propeller  stall. 

It  has  been  noted  in  the  test  data  that  the  sum  of  the 
static  pressure  in  the  inlet  and  the  dynamic  pressure  in  the 
inlet  is  not  zero.  This  indicates  that  there  is  an  inlet 
energy  loss.  To  account  for  this  loss,  an  inlet  energy  loss 
factor  was  defined  based  on  the  inlet  dynamic  pressure. 

Some  data  are  available  on  an  inlet  loss  factor  of  this  sort 
in  the  fan  literature. 

The  modification  to  the  theory  which  includes  airfoil 
data  only  partially  corrects  the  performance  for  stalling 
of  the  propeller.  When  there  is  a  large  amount  of  the 
propeller  stalling,  the  flow  breaks  down  and  can  flow  back 
through  the  propeller  plane.  The  theory  can  not  predict 
performance  for  these  conditions. 

The  modified  theory  is  presented  in  the  following 
sections. 


1.1.1  Propeller  Induced  Velocity 

The  momentum  theory  states  that  the  net  force 
acting  instantaneously  on  the  fluid  within  a  control  volume 
is  equal  to  the  time  rate  of  change  of  momentum  within  the 
control  volume  plus  the  excess  of  outgoing  momentum  flux 
over  incoming  momentum  flux.  Referring  to  Figure  1,  the 


control  volume  is  drawn  as  follows:  consider  the  zero 
streamline  (q  -  0)  to  form  one  boundary  of  the  control 
volume.  The  next  boundary  follows  the  contour  of  the 
inlet  ring  and  the  duct  walls.  Finally,  the  remaining 
boundary  is  drawn  at  the  exit  plane  of  the  GEM.  The  net 
forces  acting  on  thes  volume  are  the  propeller  and  inlet 
ring  forces,  T  and  T  ,  as  well  as  the  pressure  force  acting 


on  the  exit  plane.  The  momentum  theory  gives  the  following 
equation: 

Ts  +Tp  -  PiSj  =pVj  Sj  . 

(1) 

Similarly  for 

an  annular  area  increment, 

dTs  +  dTp  -  Pj  dS  j  *  pVj  *d  S  j . 

(2) 

Equation  (2)  can  be  rewritten  as 

/  \ 

dTp1*“  a(Pj^2qj)2‘irR/dxjxj. 

(3) 

In  order  to  express  Equation  (3)  in  terms  of  the  parameters 
behind  the  propeller,  denoted  by  the  assumption  is  made 
that  the  mass  flow  through  any  percentage  area  remains 
constant;  that  is, 


V®27TRffl*dx®x©=  Vi27rRj2dxjx 

or 

Xj dXj  _  1  -  X j  since  V©  _  ! 

x®dx© "  1-xp*  Vj  "1 

Substituting  this  into  Equation  (3) , 


ilLiiL  =/p., 

dx  1  dTp  'Kj 


Pj+2q^ 


2x 

1-XP! 


The  continuity  equation  gives 

qoSj*=q©Sr\ 

and  Bernoulli's  equation  gives 

Po+qj=P®+q©  • 


XJ 

(4) 

.  Sj  . 

(5) 

Sp 

■ 

(6) 

(7) 

(8) 
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Substituting  Equations  (7)  and  (8)  into  Equation  (6)  and 
solving  for  q®. 


since  q®  =  ipVp 


i 

z 


or  in  non  dimensional  form, 


Vp  . 

[Mg 

Spl 

[sjj 

j 

dTs  P*2x 

Si 

Sp 

Q.R  “ 

L  x 

1  + 

Sp)2 

Si 

n  1 

1  dTp  dTp  1 

1  dx  /J 

The  momentum  theorem  can  now  be  applied  to  a 
control  volume  upstream  of  the  propeller  to  give 


(10) 


(11) 


dTs  =(p,»2qi)dSi 


(12) 


or 


'dT.\ 

idTpJ 


dTp  =  (r  +2qi) 


dSi 

d  Sp 


dSp  • 


(13) 


In  order  to  account  for  inlet  losses,  it  was  found  appro¬ 
priate  to  define  an  inlet  loss  factor  related  to  the  inlet 
dynamic  pressure  as  follows: 


Pi  +  qi  -  kqi  =  0  '  (14) 

Finally,  the  propeller  thrust  is  determined,  from  a  control 
volume  enclosing  the  propeller,  as 

Tp  =(  R-P,)Sp  •  (15) 

Considering  these  relations  along  with  the  assumption  that 
anc*  sQs  =  Equation  (13)  gives  the 

following  relationship: 
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2-n-R‘xP.  ,  Ts  (Sp  (  k-1  1 

Tf  |  Si  1  k»1 
dx 

Equation  (16)  can  be  substituted  into  Equation  (11)  to 
obtain 


fl  _yDa|  dCjl  S?L  1  ,  _  I 

'  *dx(si  2  T.  Si/.,  T.rsp]fk-1l\  1 

~Trjsppj  •  <17> 

[  si)  J 

For  the  case  of  Sp=Sj=Si  ,  Equation  (17)  reduces 


v, 


An  evaluation  of  the  inlet  loss  factor  is  given  in  Section 
3.2.4. 

1.1.2  Propeller  Thrust  Coefficient 


Propeller  thrust  can  be  determined  using  the 
blade  element  theory.  Referring  to  the  following  nomem- 
clature  diagram. 


where 


dT  =  dL  cos0-dDsin0 

u^Vp^jxnR)1 

-i  {  Vp 


0  =  tQn  (xnR 

Equation  (20)  can  be  rewritten  as 


(20) 

(21) 

(22) 


and  also 


x2 

1  /Vp  l*  ~ 

z 

^  ^  / 

Vp 

(l-Xp* 

x*|nRj 41 

Ci-Cd 

xqr] 

(23) 


X  =1 

dGr 
dx 

X=  Xp 


Ct=  J ^  dx 


(24) 


1.1.3  Propeller  Power  Coefficient 


Propeller  power  can  be  written  referring  to  the 
nomenclature  diagram  as  follows: 


dP=QdQ 


(25) 


or 


dP=Qr(dDcos0+dLsinj&) 


(26) 


Similar  to  the  development  of  propeller  thrust  coefficient, 
the  power  coefficient  can  be  written  as  follows: 


dCp.  CTk  X3 
dx 


(1-Xp5 


V 


QR 


C<j+Cj 


Vp 


XQR 


(27) 
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and 


X  =1 


X=Xp  (28) 

1.2  Determination  of  the  Performance  of  a  Propeller  of 
‘  Arbitrary  Planform  and  Twist 

The  following  procedure  can  be  used  to  determine  the 
*  performance  of  an  arbitrary  propeller  at  any  operating 

condition.  Section  1.1  presents  the  propeller  power  and 
thrust  coefficient  performance  equations  along  with  the 
axial  velocity  prediction.  The  form  of  the  equations 
presented  is  such  that  they  can  best  be  solved  through  the 
use  of  numerical  integration  across  the  span  of  the 
propeller.  The  numerical  integration  involves  solving 
Equations  (17),  (23),  and  (27)  at  various  positions  on 
the  blade  span,  then  summing  the  result  to  obtain  power, 
thrust,  and  velocity.  This  process  can  be  carried  out 
in  tabular  form.  The  equations  can  be  set  up  in  a  column 
form  and  the  summation  process  can  be  carried  out.  The 
procedure  is  as  follows: 

1.2.1  Select  Radial  Stations 

Divide  the  propeller  span  into  convenient  radial 
stations  for  the  calculation  process.  The  number  of  radial 
stations  required  depends  upon  the  accuracy  desired.  A 
minimum  number  of  five  radial  stations  was  found  to  be 
adequate  to  obtain  reliable  results  for  this  program. 

1.2.2  Use  Design  Characteristics  as  Input  Data 

The  following  propeller  parameters  are  to  be 
obtained  for  use  as  input  data: 

Number  of  blades 

* 

i  Diameter 

l  Chord  (at  each  radial  station) 
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crs’-K  jasat!1"  K" 


,'*7T^'*vTSr'.^*  •iJvr  U' .  *  'V .  ■,*'  .  i.'<'Vv  'fc! 1 


Pitch  (at  each  radial  station) 

Airfoil  Section  of  blades 
Propeller  tip  speed 

A  desired  propeller  r.p.m.  is  generally  available.  If 
there  is  no  specification  for  r.p.m.,  select  a  value. 

The  lowest  r.p.m.  which  will  produce  the  desired  pressure 
without  blade  stall  generally  results  in  the  best  per¬ 
formance.  As  a  first  approximation,  a  value  of  r.p.m. 
may  be  selected  which  will  give  a  tip  speed  of  550  feet/ 
second. 


1.2.3  Assume  Axial  Velocity  Ratio,  Vp/R  R 


The  performance  equations  developed  in  Section  1 
require  an  iterative  process  for  the  solution  of  Vp/QR  at 
each  radial  station.  This  process  begins  with  the  assumption 
of  the  axial  velocity  ratio  at  the  first  radial  station 
chosen  closest  to  the  propeller  hub.  It  is  suggested  that 
a  value  of  0.15  be  chosen  to  keep  the  iteration  process  to 
a  minimum.  This  value  was  found  to  be  a  satisfactory 
starting  value  for  most  conditions. 

1.2.4  Calculate  Local  Inflow  Angle,  0 


Equation  (22) , 


The  inflow  angle  can  be  calculated  from 


0  =  tan' 


Vp 


XftR 


where  x  and  Vp/Q.  R  are  the  values  used  in  the  previous  steps. 

1.2.5  Calculate  Local  Blade  Angle  of  Attack,  OCk 

The  angle  of  attack  is  given  from  Equation  (19) 
as 


a*=e„-0 
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where  ©*  is  Che  local  blade  pitch. 


1.2.6  Calculate  Mach  Number  at  Each  Blade  Radial 
Station,  Mx 

The  Mach  number  at  each  radial  station  chosen 
from  1.2.1  can  be  calculated  as 


M 


XQR 
x"  1120 


for  standard  atmospheric  conditions. 


(29) 


1.2.7  Obtain  Section  Airfoil  Lift  and  Drag  Coeffici¬ 
ents,  Cj  and  Cd 


Consult  wind-tunnel  data  for  the  propeller  air¬ 
foil  sections.  Lift  and  drag  coefficients  can  be  obtained 
from  these  data  for  the  given  angle  of  attack  and  Mach 
number. 


1.2.8  Obtain  the  Value  of  the  Ratio  of  Inlet  Ring 

Thrust  to  Propeller  Thrust,  Ts/Tp 

The  value  of  Ts/Tp  can  be  determined  from 
Section  3.2.4. 

1.2.9  Perform  Iterative  Digital  Process 

The  information  collected  up  to  this  point  can 
now  be  used  to  calculate  performance.  The  velocity  ratio 
assumed  in  1.2.3  and  the  resulting  blade  section  lift 
coefficient  are  substituted  into  Equation  (23)  to  obtain 
dCx/dx  at  the  blade  hub  station  chosen  in  1.2.1.  This  value 
of  <jCT/dx  is  then  substituted  into  the  axial  velocity 
Equation  (17) ,  The  solution  of  Equation  (17)  represents  the 
value  of  the  calculated  axial  velocity  ratio.  This  value 
is  now  compared  with  the  assumed  value  of  1.2.3.  If  the 
values  agree,  continue  to  the  next  radial  station  and  repent 
the  process  with  1.2.3;  if  the  values  do  not  agree,  intro¬ 
duce  the  calculated  value  of  Vp/Q  R.  from  1.2.9  into  1.2.3. 
This  iterative  process  should  continue  until  the  difference 
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r 


IPIIII^HJ  II 


between  the  assumed  value  of  Vp /Q  R  and  the  calculated 
value  is  less  than  0.001.  After  satisfying  this  require¬ 
ment,  continue  to  the  next  radial  station  and  start  with 
1.2.3.  This  iterative  process  is  continued  until  complete 
radial  coverage  of  the  propeller  is  reached,  i.e.,  from 
X  -  Xp  to  X  -  1.0. 

1.2.10  Calculate  Propeller  Thrust 

Propeller  thrust  coefficient  is  obtained  by 
adding  the  incremental  thrust  coefficients,  dCx/dx,  obtain¬ 
ed  from  1.2.9  according  to  the  equation, 

X  =  1 


where 

AX  =*  non  dimensional  distance  between  each  of 
the  radial  stations  chosen. 

The  propeller  thrust  can  be  obtained  as 


X  =Xp 


dGr 

dx 


Ax 


(30) 


Tp=CtP(QR]  Sp  •  (31) 

1.2.11  Calculate  Propeller  Power 

Propeller  power  coefficient  is  obtained 
through  calculation  of  Equation  (27),  dCp/dx,  and  numer¬ 
ically  integrating  these  values  according  to  the  following 
equation: 


*£  &  - 

X  =  Xp 


(32) 
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The  propeller  shaft  horsepower  is  given  as 
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SHP  = 


Cp 


pM 


Sp 


550  (33) 

1.2.12  Calculate  Propeller  Axial  Velocity  Distribution 


The  results  of  1.2.9  give  the  exial  velocity 
ratio  at  each  radial  station.  The  velocity  distribution 
can  be  obtained  by  multiplying  the  axial  velocity  ratio  by 
the  propeller  tip  speed, Or. 


2.  DEVELOPMENT  OF  AN  OPTIMIZED  PROPELLER 


2.1  Propeller  Optimization 

2.1.1  General  Considerations 


There  are  three  general  considerations  to  be 
made  in  the  optimizing  of  a  propeller  for  a  GEM.  These 
considerations  are: 

a.  The  propeller  should  produce  the  required 
air  power  for  the  least  shaft  power  input. 

b.  The  propeller  should  produce  nearly  uniform 
velocity  at  the  jet  nozzle  of  the  GEM. 

c.  The  propeller  should  be  of  the  proper  size 
to  result  in  a  minimum  weight  of  the  pro¬ 
peller,  duct,  and  propeller  drive  system. 

This  type  of  propeller  would  lead  to  a  GEM  with  the  largest 
payload  for  a  given  mission. 

2.1.2  Specific  Considerations 

The  optimum  propeller  requirement  for  pro¬ 
duction  of  a  given  air  power  for  the  least  shaft  power 
input  can  be  expressed  as  a  propeller  efficiency  as  follows: 


V 


Air  Horsepower 
Shaft  Horsepower 


This  can  be  related  to  the  propeller  since 


T 

— P.JP.  ... 

550  IHP 


(34) 


(35) 
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or,  in  non  dimensional  form, 


77.  h  2i  • 

0R  Cp  (36) 

This  efficiency  can  also  be  expressed  as  a  propeller  blade 
section  efficiency  as  follows: 


If 


Vp  j  dCr 

7?*-- ic—  ■ 

dx 


it  is  assumed  that 


Cd(Vp/xnRj«  C, 

{vP/  XQR  ) 2  <(  1.0 


and 


Sp  =Sj  =Si 


(37) 


then  substituting  Equations  (18),  (23),  and  (27)  of 
Section  1.1  into  Equation  (37)  of  this  section,  it  is  found 
that 


_ 1 

1 

'  crx  I _ t; 

[2x(Tp|U 
This  equation  indicates  that  i: 
will  be  defined  as  the  airfoil  section  efficiency  parameter, 
is  increased,  the  propeller  section  efficiency  will  increase. 
This  increase  occurs  in  a  fairly  complicated  manner,  and 
therefore  the  following  typical  values  were  assumed  to  evalu* 
ate  this  relation: 


,  .£*■ 

T  Cd 


•X 

r 


(38) 


the  ratio  Cj3/2/Cd,  which 
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0.22 


o; 


x  -  0.80 
Ts/Tp  -  0.20 
k  -  -0.53 

The  blade  section  efficiency,  calculated  for 
various  values  of  the  airfoil  section  efficiency  parameter, 
is  shown  in  Figure  2.  It  should  be  noted  from  this  figure 
that  a  small  increase  in  efficiency  results  from  a  large 

increase  in  Cj'^/Cd  when  this  parameter  is  above  a  value 
of  40. 


It  may  also  be  noted  from  the  relation  for 
blade  section  efficiency  that  the  blade  section  solidity, 
Ox  and  Ta/Tp,  should  be  increased  for  an  increase  in 

efficiency.  Also,  the  inlet  loss  factor,  k  should  be 
small. 


The  optimum  propeller  requirement  for  pro¬ 
duction  of  uniform  velocity  to  produce  best  GEM  performance 
was  not  analyzed  further  in  this  program.  However,  it  should 
be  noted  that  the  analysis  presented  in  Reference  1  showed 
that  GEM  performance  was  not  significantly  reduced  unless 
the  velocity  differed  greatly  from  that  which  is  uniform. 

This  result  was  substantiated  by  the  test  program  as  noted 
later  in  this  report.  As  a  consequence  it  may  be  found 
that,  in  the  evaluation  of  an  optimum  propeller,  other  con¬ 
siderations,  such  as  weight  or  cost,  may  be  more  significant 
than  the  effect  of  producing  a  uniform  velocity.  These 
factors  were  not  considered  in  this  program. 

To  satisfy  the  requirement  for  a  propeller, 
duct,  and  drive  system  of  minimum  weight  requires  a  design 
study  of  the  specific  configuration.  Therefore,  it  is  not 
possible  to  generalize  as  to  the  size  of  the  propeller. 
However,  the  requirement  for  compact  ducting  with  small 
diffusion  loss  indicates  that  a  propeller  area  of  about  the 
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same  size  as  the  GEM  nozzle  area  is  required.  Also,  the 
propeller  size  which  results  in  an  adequate  propeller  tip 
speed  without  a  transmission  between  the  engine  and  propeller 
is  attractive.  On  the  other  hand,  large-area  propellers 
tend  to  occupy  the  usable  space  on  the  GEM.  These  consider¬ 
ations  are  beyond  the  scope  of  the  present  program. 

2.2  Optimized  Propeller  Design  Procedure  and  Considerations 

The  design  procedure  for  an  optimized  propeller  is  out¬ 
lined  in  this  section,  along  with  information  which  is 
believed  to  be  pertinent  to  the  design  considerations.  This 
design  procedure  is  enumerated  in  detail  with  the  design  of 
the  optimized  propeller  used  in  the  test  program  as  an 
example  in  Appendix  III. 

It  is  to  be  noted  that  some  of  the  procedures  presented 
in  this  section  are  interchangeable.  That  is,  should  one 
part  of  the  design  of  the  propeller  characteristics  be 
already  established  (for  example,  the  propeller  diameter), 
it  is  possible  to  incorporate  this  design  feature  with  the 
remaining  design  steps. 

The  performance  equations  developed  in  Section  1.1  are 
used  for  determining  the  characteristics  of  the  optimized 
propeller.  Thus,  the  procedure  is  similar  to  the  develop¬ 
ment  of  Section  1.2,  which  presents  the  method  of  determin¬ 
ing  arbitrary  propeller  performance.  As  was  mentioned  in 
Section  1.2,  the  form  of  the  equations  presented  is  such 
that  they  can  best  be  solved  through  the  use  of  a  numerical 
integration  across  the  span  of  the  propeller. 

The  first  part  of  the  design  procedure  involves 
obtaining  an  optimum  blade  airfoil  section.  An  optimum 
airfoil  as  described  in  Section  2.1  is  one  which  produces 

maximum  for  the  expected  operating  conditions. 

Thus,  a  study  of  airfoil  data  is  required  to  reveal  such  an 
airfoil.  The  angle  of  attack  corresponding  to  C*3/2/Cd 
maximum  can  then  be  obtained  for  the  airfoil  selected 
at  various  Mach  numbers. 


17 


Propeller  tip  speed, &R,  is  then  chosen  to  insure 
operation  compatible  with  Mach  number  and  structural  consid¬ 
erations.  Propeller  tip  speed  variation  up  to  600  f.p.s. 
was  noted  to  have  little  effect  upon  performance  and  the 
optimization  process.  The  lower  limit  to  tip  speed  results 
from  structural  considerations  and  from  the  large  coning 
angle  which  must  be  built  into  the  propeller  as  explained 
in  Appendix  II.  The  high  limit  to  practical  tip  speeds  is 
caused  by  compressibility  effects.  For  thin  airfoil  sect¬ 
ions  at  low  angles  of  attack,  satisfactory  performance  can 
be  achieved  up  to  about  900  f.p.s.  However,  the  range  of 
tip  speeds  from  400  to  600  f.p.s.  is  recommended  for  initial 
design  purposes. 

The  following  area  parameters  are  to  be  obtained: 

a.  Propeller  to  inlet  area  ratio,  Sp/Si 

The  effect  of  this  parameter  upon  the  optimiz¬ 
ation  process  was  not  studied  during  this  program. 

It  is  believed,  however,  that  an  inlet  area  smaller 
than  the  propeller  area  would  reduce  inlet  losses, 
since  the  smaller  inlet  area  would  tend  to  increase 
inlet  dynamic  pressure,  thus  preventing  separation. 

b.  Propeller  to  jet  area  ratio,  Sp  /  Si 

Tests  indicated  that  the  best  GEM  performance 
was  obtained  for  the  case  in  which  the  propeller 
area  was  equal  to  the  jet  area  for  the  range  of  h/d 
from  0.05  to  0.30.  Also,  the  smallest  inlet  losses 
occurred  for  the  case  in  which  Sp/Sj  "1.0 

Propeller  axial  velocity  can  be  obtained  through  the  use 
of  GEM  design  parameters  as  follows: 

a.  An  intial  design  study  will  usually  lead  to  such 
parameters  as  vehicle  gross  weight,  size,  and 
maximum  operating  height  required. 
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b.  The  average  base  pressure  required  to  support  the 
given  weight  for  this  vehicle  is 


Pb=  W 

rb  s 


(39) 


c.  GEM  nozzle  geometry,  such  as  jet  area  and  inclina¬ 
tion  angle,  is  to  be  obtained.  Reference  2  presents 
the  effects  of  jet  turning  angle  upon  performance. 

d.  The  ratio  of  base  pressure  to  jet  total  pressure, 
Pb/Ptj »  and  the  ratio  of  the  static  pressure  to  the 

total  pressure  at  the  jet,  Pj/Pt  >  can  be  obtained 

for  the  given  nozzle  geometry  and  operating  height 
from  various  sources.  In  Reference  1,  a  compilation 
of  some  of  the  theoretical  analysis  along  with  test 
data  is  presented. 

e.  The  jet  velocity  can  be  calculated  from  the  given 
ratios  as 


qj^Ptj-Pj  • 


(40) 


This  equation  can  be  reduced  in  terms  of  the  data 
previously  obtained  as 


and  also 


q  .  =  _pb 

HJ  _Pb 
Ptj 


Pj 

Ptj) 

Iptj 

Pb 

(41) 


V4  =■/ 


2. 

P 


qj 


(42) 


f.  Propeller  axial  velocity  can  be  calculated  from 
the  continuity  equation  as 
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(43) 


Vp=  Vj 


The  condition  of  constant  propeller  axial  velocity 
ratio,  Vp/Q R,  can  be  imposed  as  a  condition  to 

solve  the  performance  equations. 

The  local  blade  pitch  as  a  function  of  blade  radius 
(that  is,  the  blade  twist)  can  be  determined  by  adding  the 
optimum  airfoil  section  angle  of  attack  to  the  local  inflow 
angle. 

Thus 


where 


Op  =  |C()  optimum  *0 


(44) 


0  =  tan'1 


VP 


XQR 


The  local  propeller  solidity  can  be  obtained  through 
the  elimination  of  dCT/dx  in  Equation  (17)  through  the  use 
of  Equation  (23)  to  give 


where  the  values  of  Ts/Tp  and  K  are  discussed  in  Section  3.2.2. 


The  solidity,  Ox  ,  varies  as  the  product  of  the  chord 
and  number  of  bl.ades  for  a  fixed  propeller  diameter.  Thus, 
when  establishing  the  number  of  blades,  consideration  must 
be  given  to  satisfying  the  solidity  with  practical  values  of 
the  chord.  This  requirement  is  generally  a  structural  problem. 
Aerodynamic  efficiency  usually  will  be  slightly  improved  with 
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a  larger  number  of  blades.  However,  the  advantages  of 
having  a  minimum  of  blades  are  brought  out  in  the  simpli¬ 
city  of  both  hub  and  control  design. 

The  number  of  blades  is  related  to  the  local  solidity 
and  the  local  chord  of  the  blades  by  the  following  relation: 


C*  = 


Cy'TT  R 
b 


(46) 


If  it  is  found  that  the  local  chords  of  the  propeller  are 
such  that  the  relation  zrrr  is  less  than  the  value  2, 

be  ’ 

cascade  effects  will  be  significant.  These  cascade  effects 
influence  the  aerodynamic  performance  of  the  blade  sections 
and  may  result  in  a  poor  prediction  of  performance.  Cascade 
correction  factors  are  presented  in  Reference  3,  which  appears 
to  be  applicable  for  these  conditions. 


3.  TEST  RESULTS 


In  this  section  the  results  of  the  test  program  are 
presented.  The  analysis  has  been  utilized  to  isolate  the 
performance  components;  however,  the  comparison  of  test 
results  with  predicted  performance  is  presented  in  Section  4. 

3.1  Propeller  Performance 


The  thrust  and  power  coefficients  of  a  propeller  define 
the  performance  of  that  propeller  for  given  operating  condi¬ 
tions.  The  thrust  coefficient  is  a  measure  of  the  increase 
in  pressure  which  the  propeller  causes.  Since  the  air  (or 
useful  power)  caused  by  the  propeller  is  a  function  of  the 
3/2  power  of  the  thrust  coefficient,  an  expression  for  the 
propeller  efficiency  can  be  obtained  from  the  performance 
coefficients.  The  propeller  data  obtained  in  this  program 
are  presented  in  this  section.  A  comparison  of  these  data 
with  the  analysis  of  Section  1  is  presented  in  Section  4. 


follows: 


3.1.1  Propeller  Efficiency 

The  propeller  efficiency  can  be  expressed  as 
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This  relation  can  also  be  expressed  as 
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Propeller  efficiencies  calculated  with  this 
relation  are  given  in  Figure  3  for  the  straight  bladed 
propeller  and  for  the  optimized  propeller  at  two  pitch 
settings.  Maximum  efficiencies  on  the  order  of  80-percent 
are  reached  atTTdh  of  1.75.  It  can  be  seen  from  this 

figure  that  the  optimized  propeller  gave  better  performance 
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at  low  values  of 
propeller  showed 


TTdh 
Sj 

better  values  at  ~ 


however,  the  straight- b laded 

larger  than  0.70. 


The  performance  of  the  optimized  propeller  with  an 
increased  pitch  of  2  degrees  (©0“  19.4  degrees)  shows 
improvement  at  large  values  of  TTdh  ,  at  the  expense  of 
poorer  performance  at  lowTTdh  .  Sj 

Sj 


3.1.2  Thrust  Coefficient 


Propeller  thrust  coefficient  data  for  the 

various  jet  area  configurations  are  shown  in  Figure  4. 

Sj 

The  configuration  with  s  of  0.8  is  seen  to  show 
the  lowest  Ct  for  the  heights  tested;  also,  the  thrust 
coefficient  is  sensitive  to  small  changes  in  height,  while 

with  Si/s  of  0.3,  Ct  is  almost  independent  of  -?r^— - 
J  Sj 

3.1.3  Power  Coefficient 


Propeller  power  coefficient  test  data  are  shown 
in  Figure  5  for  the  various  jet  area  configurations  tested. 
It  may  be  noted  that  the  variation  of  this  parameter  with 
height  is  similar  to  the  thrust  coefficient  data  in  that 
power  coefficient  is  independent  of  height  for  the  small 
jet  area  configurations. 

3.2  Duct  Performance 


The  GEM  duct  system  can  be  divided  into  two  parts  for 
the  purpose  of  analysis.  The  first  portion  of  the  system 
in  front  of  the  propeller  consists  of  an  inlet  ring  and  an 
inlet  duct;  the  second  part  of  the  system  contains  the  duct 
from  the  propeller  to  the  exit  nozzles.  In  this  section,  the 
performance  of  the  complete  duct  system  is  considered  and 
then  the  inlet  performance  is  discussed  separately. 

3.2.1  Duct  Efficiency  and  Loss  Factor 

The  efficiency  of  the  duct  system  of  the  GEM  test 
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unit  is  defined  to  be  the  ratio  of  nozzle  air  power  to  the 
air  power  produced  by  the  propeller, 


7)  -  Vj 

/o  '  TpVp 

or,  by  use  of  the  continuity  equation,  Si  Vi  “SpVpj 


(49) 


(50) 


The  duct  efficiency  was  calculated  with  the 
measurement  of  total  jet  pressure,  Ptj ,  and  propeller  thrust, 

Tp,  according  to  Equation  (50)  Figure  6  gives  Vo  for  various 
configurations  for  the  optimized  and  the  straight-bladed 
propeller.  The  average  efficiency  of  76- p-  .  ,ent  is  seen  to 
be  practically  independent  of  conf igurati^u  and  nozzle  height 
for  the  optimized  propeller.  The  duct  efficiency  with  the 
straight-bladed  propeller  is  4  to  6-percent  higher  than  the 
duct  efficiency  with  the  optimized  propeller.  This  increased 
efficiency  is  caused  by  the  increase  in  the  inlet  efficiency. 
The  straight-bladed  propeller  produced  a  velocity  distribution 
which  increased  toward  the  blade  tip.  This  increased  velocity 
near  the  duct  wall  prevented  inlet  ring  flow  separation,  thus 
increasing  inlet  efficiency  and,  hence,  duct  efficiency. 


The  efficiency  of  a  duct  system  downstream  of 
the  propeller  can  also  be  defined  in  terms  of  dynamic  pres¬ 
sure  losses.  A  duct  loss  factor  kz  is  defined  by  the  follow¬ 
ing  relation: 


P.*q.(l-k!)  =  Pj.qj  •  (5! 

The  duct  efficiency  of  Equation  &9) can  be  expressed  in  terms 
of  the  GEM  nozzle  static  total  pressure  ratio  with  the  use 
of  Equation  (51) and  the  continuity  equation  as  follows: 
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The  duct  loss  factor  is  shown  to  approach  infinity 
as  sfj — approaches  zero.  This  is  due  to  the  fact  that  the 

jet  static  pressure  becomes  equal  to  the  jet  total  pressure 
at  zero  height.  The  large  values  of  the  duct  loss  factor 
increasing  with  decreasing  height  and  decreasing  jet  area  are 
believed  to  be  caused  by  the  inlet  losses. 

3.2.2  Inlet  Thrust  and  Inlet  Loss  Factor 

Inlet  thrust  is  produced  by  the  negative  static 
pressure  acting  on  the  inlet  ring  area.  The  magnitude  of  the 
inlet  static  pressure  is  dependent  upon  GEM  height.  This  in¬ 
let  thrust  was  determined  through  the  use  of  measured  surface 
static  pressures  and  by  use  of  strain  gages  as  described  in 
Appendix  I. 

The  ratio  of  inlet  ring  thrust  to  propeller  thrust 
for  the  straight-bladed  and  the  optimized  propeller  is  shown 
in  Figure  8  as  a  function  of  1T.d._h_.  .  The  two  propellers 

Ts  Sj  TTd  h 

show  a  similar  -  variation  for  values  of-1— —  below  3.0. 

Ip 

The  parameter  Sj/S  is  seen  to  have  a  pronounced  effect  upon 

Jk  • 

Tp 

As  discussed  in  Section  1.1.1,  an  inlet  loss  factor 
can  be  defined  as 


Pi  +  qi  =  kq,  •  (53) 
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The  relation  between  inlet  ring  thrust  and  propeller  thrust 
can  thus  oe  written  in  terms  of  duct  pressures  and  the 
inlet  loss  factor  as  follows: 
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The  toss  factor,  k,  is  available  from  the  literature  in  some 
cases.  Values  of  k  were  determined  with  the  use  of  Equation 
(54)  and  from  the  test  data.  The  inlet  loss  factor  would 
be  equal  to  zero  for  a  100-percent  efficient  inlet.  Figure 
9  presents  the  inlet  loss  factor  for  various  jet-to-base- 
area  ratios.  It  is  seen  from  this  figure  that  the  factor 
approaches  zero  with  increasing  values  of  .  The  larger 

jet-area  configurations  thus  increase  the  inlet  efficiency 
by  increasing  the  GEM  volume  flow  at  the  same  TT  d  h  . 

3  j 

3.3  GEM  Performance 


One  of  the  main  topics  of  study  in  this  program  was 
the  determination  of  the  effect  of  the  optimized  propeller 
on  GEM  performance.  It  was  noted  from  the  analysis  of 
Reference  1  that  the  uniform  nozzle  velocity  caused  by  the 
optimum  propeller  would  result  in  improved  performance  as 
compared  with  the  straight  Dladed  propeller.  The  GEM  per¬ 
formance  test  data  obtained  with  these  two  propellers  are 
compared  in  the  following  discussion. 

The  effect  on  GEM  performance  of  varying  the  ratio  of 
the  GEM  nozzle  (jet)  area  to  the  base  area  is  also  discussed 
in  this  section. 

It  should  be  noted  that  GEM  performance  (that  is,  lift 
and  nozzle  air  power  required)  can  be  determined  if  Pb/Pti 
and  Pj/Ptj  are  known.  J 

3.3.1  Base  Pressure  -  Nozzle  Total  Pressure  Ratio 

The  data  obtained  on  the  base  pressure  nozzle 
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total  pressure  ratio  are  shown  in  Figure  10.  It  may  be 
noted  in  these  data  that  there  is  no  effect  of  the  pro^ 
peller  configurations  on  these  data.  The  effect  of 
on  this  ratio  at  a  given  value  of  ■■■■—■ J*. - ■■■  is  small, 
but  a  slight  increase  in  may  Jbe  noted  for 

Sj/S  of  0.8.  P*j 

For  comparison  purposes,  the  test  data  of 
Reference  4  are  also  shown  on  Figure  10.  The  GEM  tested 
in  this  reference  had  an  Sj/S  of  0.13  and  a  design  nozzle 
inclination  of  30  degrees.  These  data  show  a  higher  base 
pressure  to  nozzle  pressure  ratio  than  the  configurations 
tested  in  this  program  due  to  the  inclination  of  the  nozzles. 

3.3.2  Nozzle  Static  -  Total  Pressure  Ratio 


The  ratio  of  jet  static  pressure  to  jet  total 
pressure  is  shown  in  Figure  11  as  a  function  of  — — 

The  curves  show  only  slight  variation  with  Sj/S.  J 

3 . 4  Overall  Performance  of  Propeller  and  GEM 

The  combined  performance  of  the  GEM  with  the  propeller 
and  duct  system  is  discussed  in  this  section.  Performance 
data  are  presented  in  dimensional  form  as  measured  on  the 
test  unit.  The  internal  efficiency  of  the  unit,  including 
propeller  and  duct  efficiencies,  is  discussed.  Overall 
performance  of  the  unit  as  a  GEM  is  presented  in  the  form 
of  the  nondimensional  power  factor,  2  . 

3.4.1  Dimensional  Performance  Data 

The  performance  of  the  GEM  test  unit  is  given  as 
the  total  lift  divided  by  input  shaft  horsepower,  L/lHP,for 
various  heights  in  Figure  12.  Power  loadings  are  shown  for 
the  straight-bladed  propeller  and  for  the  optimized  pro¬ 
peller  with  various  nozzle  configurations.  The  optimized 
propeller  test  data  show  a  considerable  increase  in  the 
pounds  per  horsepower  for  similar  GEM  lift- total  area  load¬ 
ings,  L/S,  for  a  range  of  h/d  from  0.010  to  0.10  over  the 
straight-bladed  propeller. 
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A  power  loading  of  21.5  pounds/horsepower  pro¬ 
duced  by  the  optimized  propeller  at  an  h/d  of  0.02  is  seen 
to  be  almost  double  that  of  the  straight-bladed  propeller. 

The  improvement  in  performance  due  to  the  change  in  the 
propellers  is  due  to  an  increased  resistance  to  stalling 
of  the  optimized  propeller. 

Figure  12  also  shows  the  effect  of  the  nozzle 
area  to  total  base  area  parameter,  Sj/S,  upon  performance 
for  the  optimized  propeller.  Similar  performance  for 
Sj/S  of  0.8  and  0.6  is  noted  for  the  range  of  h/d  tested, 
while  lower  power  loadings  are  noted  for  Sj/S  of  0.3  and  0.10. 

3.4.2  Internal  Efficiency 

The  internal  efficiency,  defined  as  the  ratio 
of  nozzle  air  horsepower  to  input  shaft  horsepower  expressed 
as  a  percentage,  includes  the  efficiencies  of  the  propeller 
and  the  duct.  Figure  13  presents  the  internal  efficiency 
forthe  straight-bladed  propeller  as  well  as  for  the  optimized 
propeller  for  various  duct  configurations.  The  straight- 
bladed  propeller  shows  higher  internal  efficiencies  for  h/d 
greater  than  0.10,  while  the  optimized  propeller  indicates 
improved  performance  for  h/d  less  than  0.10.  The  internal 
efficiency  decreases  with  decreasing  jet  area;  however,  it  is 
noted  that  for  the  case  of  Sj/S  equal  to  0.6,  the  internal 
efficiency  is  higher  for  h/d  greater  than  0.15  than  for  the 
data  obtained  with  Sj/S  of  0.8.  This  slight  improvement  is 
due  to  the  increased  propeller  efficiency  at  h/d  greater  than 
0.15. 

3.4.3  GEM  Power  Factor,  Z 


In  order  to  evaluate  the  effect  of  the  propeller 
combination  and  the  jet  area  relation  on  GEM  performance,  it 
is  appropriate  to  define  an  overall  power  factor  as  the  ratio 
of  input  shaft  horsepower  to  the  horsepower  required  by  an 
ideal  ducted  fan  hovering  out  of  ground  effect,  with  a  disc 
loading  equal  to  W/S,  or  in  equation  form, 
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This  equation  can  also  be  expressed  in  terms  o£  the  measured 
internal  efficiencies  of  Section  3.4.2  and  the  calculated 
volume  flows  as 


2PtJF 
Vi  L  i 

ips 
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Figure  14  presents  fKis  factor  as  a  function  of  the  nozzle 
height  parameter  £o:  various  values  of  Sj/S.  It  should  be 
noted  that  the  GEM  performance  increases  with  decreasing 
values  of  the  overall  power  factor.  The  optimized  propeller 
shows  improved  performance  over  the  straight- bladed  propeller 
for  values  of  h/d  less  than  0.3.  It  is  also  noted  that  the 
jet-base  area  ratio  of  0.6  configuration  shows  increased 
performance  over  jet-base  area  ratio  of  0.8  configuration 
for  h/d  values  greater  than  0.10. 
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4.  COMPARISON  OF  PREDICTED  PERFORMANCE  WITH  TEST  RESULTS 


In  this  section,  the  performance  predicted  by  the  analy¬ 
sis  of  Section  1  is  compared  to  the  test  data  presented  in 
Section  3.  An  analysis  of  these  results  is  made  to  determine 
why  the  theory  does  not  give  a  better  prediction  of  perform¬ 
ance. 


4.1  Thrust  and  Power  Coefficients 

~’*~'Hcted  and  experimental  thrust  and  power  co- 

effic- _  v.  compared  in  Figures  15,  16,  17,  and  18 

for  the  four  Sj.  _  configurations  tested.  These  data  show 
that  the  inlet  loss  theory  of  Section  1  gives  a  good  pre¬ 
diction  of  the  thrust  coefficient  and  power  coefficient  for 
Sj/S  of  0.8.  The  predicted  power  coefficient  differs  from 
the  experimental  values  at  low  height;  however,  the  pre¬ 
diction  is  within  12- percent  of  the  experimental  values, 
even  at  'Ir —  of  0.4.  As  shown  in  Figure  15,  the  inlet 
loss  theory  from  Section  1  predicts  the  experimental  data 
with  more  accuracy  than  the  theory  of  Reference  1.  The 
predicted  power  and  thrust  coefficient  for  the  smaller  jet 
area  configurations  are  not  in  agreement  with  the  experi¬ 
mental  data  as  seen  from  Figures  16,  17,  and  18. 

The  difference  between  experimental  data  and  predicted 
values  is  due  to  the  poor  prediction  of  the  propeller  blade 
section  angle  of  attack  and  blade  stall  effects,  as  will  be 
discussed  in  the  next  section. 

4.2  Propeller  Section  Angle  of  Attack 

An  estimate  of  the  angle  of  attack  of  the  propeller 
blades  can  be  made  by  assuming  that  the  velocity  at  the 
propeller  is  uniform.  With  this  assumption,  the  nozzle 
flow  data  from  Figures  6  and  11  were  utilized  to  prepare 
estimated  angles  of  attack  at  the  root  of  the  optimized 
propeller  for  various  configurations.  These  data  are  com¬ 
pared  with  the  angles  of  attack  predicted  by  the  inlet  loss 
theory  of  Section  1  in  Figure  19.  The  root  angles  of  attack 
are  shown  since  root  stall  occurs  before  blade  tip  stall, 


30 


X. 


due  to  the  blade  twist  of  the  optimized  propeller.  Blade 
stall  angle  of  attack  which  was  determined  from  the  data 
of  Reference  6  occurs  at  the  angle  of  attack  at  which  the 
lift  curve  fails  to  increase  with  an  increase  in  angle  of 
attack.  This  angle  was  found  to  be  9  degrees.  It  may  be 
noted  in  Figure  19  that  the  root  of  the  propeller  blade 
was  stalled  for  Sj/S  of  0.3  and  0.1  for  all  of  the  heights 
tested.  The  theory  predicted  that  there  would  not  be  stall 
for  these  conditions.  The  accuracy  of  the  prediction  of 
blade  section  angles  of  attack  has  been  studied  in  more 
detail,  and  the  results  are  shown  in  Figures  20  and  21.  The 
radial  distribution  of  angle  of  attack  for  Sj/S  of  0.8  is 
shown  in  Figure  20  for  two  heights.  A  very  good  prediction 
of  angle  of  attack  is  indicated  by  these  data,  especially 
for  the  inlet  loss  theory.  In  Figure  21,  similar  data  are 
shown  for  the  Sj/S  of  0.6.  It  may  be  noted  that  due  to 
stalling,  very  large  angles  of  attack  were  measured  at  the 
root  section  of  the  blade  at  a  i  h —  of  0.34.  The  theory 

predicted  small  angles  of  attack  and  no  blade  stall  for 
this  condition. 
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APPENDIX  I 


THE  TEST  PROGRAM 


1.0  SCOPE  AND  PURPOSE 

The  purpose  of  this  program  was  to  determine  the 
extent  of  applicability  of  the  analysis  of  Reference  1 
in  the  determination  of  an  optimum  propeller  for  Ground 
Effect  Machines.  A  propeller  (herein  referred  to  as  the 
optimized  propeller)  was  designed  using  this  analysis  and 
was  tested  at  various  GEM  nozzle- to-ground  heights  and  at 
two  propeller  disc  loadings.  Additional  tests  were  also 
made  with  the  straight-bladed  propeller  tested  in  Reference 
1  to  insure  that  the  data  obtained  were  comparable.  The 
optimized  propeller  was  tested  at  two  pitch  settings.  The 
effect  of  the  ratio  of  the  base  area  to  jet  area  upon  the 
optimization  process  was  also  investigated. 

2.0  DESCRIPTION  OF  TEST  APPARATUS 

The  GEM  test  device,  shown  in  Figures  22  and  23,  con¬ 
sisted  of  a  propeller  and  drive  system,  an  inlet  ring,  a 
duct,  and  a  ground  board,  instrumentation  was  provided  to 
measure  GEM  performance  and  flow  parameters. 


2.1  Propellers 


The  characteristics  of  the  propellers  tested  in 
this  program  are  shown  in  Figures  24  and  25.  The  straight 
bladed  propeller  had  a  constant  chord  and  a  0012  airfoil 
section.  The  optimized  propeller  shown  in  Figure  25  was 
twisted  and  tapered  and  had  a  16-709  airfoil  section. 

Both  propellers  had  four  blades,  were  of  wooden  construc¬ 
tion,  and  had  a  solidity  of  0.22.  As  noted  in  Figure  25, 
the  optimized  propeller  was  tested  at  root  pitch  settings 
of  17.5  and  19.4  degrees. 


2.2  Drive  System 


Power  was  delivered  to  the  propeller  by  a  hydrau¬ 
lic  drive  system  which  was  driven  by  a  20-horsepower  elect¬ 
ric  motor.  The  hydraulic  drive  permitted  a  wide  range  of 
adjustment  of  the  propeller  speed  to  produce  the  required 
propeller  disc  loadings. 

Propeller  thrust,  r.p.m.,  and  torque  instrumenta¬ 
tion  were  provided  in  the  drive  system. 

2.3  Inlet  Ring 

The  inlet  ring  had  a  semicircular  cross  section 
with  a  2.25-inch  radius,  and  an  inside  diameter  of  2.90  feet 
to  match  the  diameter  of  the  duct  svstem.  This  inlet  is 
shown  on  the  test  device  in  Figures  22  and  23. 

2.4  Duct  System 

The  duct  consisted  of  a  circular  straight  duct 
2.90  feet  in  diameter  by  4.50  feet  long,  mounted  concentric 
to  various  centerbodies  which  housed  the  drive  unit.  The 
centerbody  was  supported  by  four  struts  which  were  of  20- 
percent  thickness  and  13- inch  chord,  as  shown  in  Figure  23. 
The  exit  nozzle  area  of  the  system  was  varied  by  changing 
the  centerbody.  Figure  23  shows  the  duct  configuration 
for  the  four  centerbodies  tested.  The  parameters  of  the 
duct  system  are  listed  in  Table  1. 

2.5  Ground  Board 


A  9-foot-diameter  disc  was  used  to  simulate  the 
ground.  This  disc  was  suspended  normal  to  the  duct  axis 
on  three  cantilever  beam  assemblies,  permitting  the 
measurement  of  differential  axial  force.  Variation  of 
the  angles  of  pitch  and  roll  as  well  as  the  distance  from 
the  duct  exit  plane  was  accomplished  by  means  of  screw 
adjustments.  Tests  were  performed  at  duct  exit  heights 
ranging  from  0.062  inch  to  10  feet  and  with  roll  and  pitch 
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angles  of  t  10  degrees. 


2.6  Instrumentation 

2.6.1  Pressure  Data 


Pressure  data  were  obtained  at  various 
stations  throughout  the  test  unit.  Pressure  probes  and 
rakes  sensed  total  and  static  pressure.  The  probes  were 
connected  by  plastic  tubing  to  four  banks  of  multiple  tube 
manometers.  Pressure  data  were  recorded  at  the  following 
locations: 


2. 6. 1.1  Inlet  Ring 

The  axial  distribution  of  static 
pressure  along  the  inlet  surface  was  measured  at  four 
azimuth  positions  of  the  inlet  ring. 

2 . 6 . 1 . 2  Duct  Pressures 

Static  and  total  pressures  were 
usually  measured  at  four  azimuthal  and  nine  radial  positions 
and  at  one  station  before  the  propeller  and  at  two  stations 
downstream  of  the  propeller  plane.  Duct  surface  static 
pressures  were  measured  along  the  duct  wall  and  along  the 
centerbody.  During  certain  tests  a  calibrated  pressure 
probe  was  used  to  obtain  data  concerning  the  propeller 
swirl  angle. 


2. 6. 1.3  Base  Pressures 


The  radial  variation  of  static 
pressure  across  one  diameter  of  the  bases  of  the  center- 
bodies  was  measured. 


2. 6. 1.4  Pressure  Surveys 

To  supplement  the  duct  pressure 
probes,  pitot  tube  surveys  were  made  in  certain  areas.  In 
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particular,  the  flow  pattern  near  the  inlet  and  also  near 
the  exit  nozzle  was  studied  at  various  azimuthal  and  radial 
positions. 


2.6.2  Performance  Data 


In  addition  to  pressure  data,  the  follow¬ 
ing  performance  data  were  obtained: 

2. 6. 2.1  Propeller  Thrust 

The  propeller  thrust  was  measured 
with  strain  gages  on  a  tension  link  attachment  of  the  drive 
shaft.  A  four-arm,  temperature-compensated  strain-gage 
bridge  circuit  was  utilized.  The  signal  was  fed  through  a 
bridge  balance  and  calibrating  unit,  and  subsequently  was 
measured  by  a  millivoltmeter. 

2. 6. 2. 2  Propeller  Torque 

The  propeller  torque  was  measured 
using  strain  gages  on  a  cantilever  beam  system  mounted  in 
the  drive  shaft.  A  circuit  similar  to  that  of  the  propeller 
thrust  measuring  device  was  used  with  the  signal  being 
measured  by  a  microammeter. 

The  propeller  torque  was  also 
obtained  by  using  a  calibrated  hydraulic  pressure  gage 
mounted  on  the  drive  unit. 

2.6.2. 3  Propeller  Rotational  Speed 

A  magnetic  pickup  which  sensed  the 
passage  of  a  propeller-shaft-mounted  gear  was  used  in  re¬ 
cording  the  propeller  r.p.m.  A  stroboscope  was  also  used  to 
obtain  a  second  measurement  of  r.p.m. 

2. 6. 2. 4  Total  GEM  Lift 

The  ground  board  was  instrumented 
on  the  three  flexures,  mentioned  previously,  using  strain- 
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gage  bridge  circuits  similar  to  those  used  for  propeller 
torque  and  thrust  measurements.  From  the  three  reaction 
loads  obtained,  the  total  GEM  lift  and  the  center  of 
pressure  of  lift  were  determined. 

2.7  Evaluation  of  Test  Results 

In  order  to  evaluate  the  performance  of  the  GEM 
and  propeller,  certain  parameters  are  needed.  The  follow¬ 
ing  sections  present  methods  of  calculating  some  of  the 
important  parameters. 

2.7.1  Performance  Data 

The  following  data  were  obtained  directly 
from  recorded  data: 


Propeller  thrust,  Tp,  pounds 

Total  lift,  L,  pounds  (Total  ground  board 
reaction) 

Propeller  speed,  N,  r.p.m. 

Propeller  torque,  Q,  foot-pounds 

Utilizing  the  recorded  data,  the  shaft 
horsepower  was  calculated  as 


IHP  = 


5250 


Also,  the  power  loading,  L/lHP,  was  calculated, 


2.7.2  Pressure  Data 


To  present  the  pressure  data  in  a  more 
usable  form,  these  data  were  evaluated  by  using  the  follow¬ 
ing  relations: 


2. 7. 2.1  Internal  Efficient 


The  internal  efficiency  is  defined 
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T?  T7 - 


as  the  ratio  between  the  air  horsepower  at  the  duct  exit 
AHP,  and  the  propeller  shaft  power  inpuc,  IHP.  That  is, 

77.  =  ah  p  , 
h  IHP 

2. 7. 2. 2  Duct  Exit  Air  Power 


The  duct  exit  (or  jet)  air  horse 
power  is  defined  by  the  following  equation: 

ahp  =  35T-  ( 

Since  Fj  and  Ptj  vary  locally  over  the  duct  exit,  this 

calculation  was  accomplished  by  a  numerical  integration 
of  the  product  of  the  area  total  pressure  and  velocity 
over  the  duct  annulus  at  the  nozzle.  Total  pressure  was 
obtained  directly  from  the  pressure  data  at  the  nozzle. 

The  local  velocities  were  calculated  from  the  static 
and  total  pressures  using  the  Bernoulli  equation. 


qj=Ptj-Pj  =  -|pvj 


Solving  for  velocity, 


Vi=  Vp-(P  tJ-PJ 


The  resulting  equation  to  be  solved  by  numerical  integra¬ 
tion  is  X-i 

ahp=Ho/V  I- (ptj'pj)ptj xdx  • 

X=XH 

The  air  density,  P,  was  corrected  for  duct  temperature 
before  use  in  this  relation. 
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2. 7. 2. 3  Volume  Flow.  F 


Volume  flow  was  calculated  by  per¬ 
forming  a  numerical  integration  of  the  velocity  over  the 
annular  duct  area, 


Fj  =VjAj 


Using  Equation  (3) ,  the  integral  form  of  Equation  (7) 
becomes  _  . 

A  *•  I 

Fs27rRzJiJ  -p-(ptj-Pj  Xdx  • 

X  =  Xm 

2.1 .2  A  Jet  Momentum, 


culated  by 


The  momentum  of  the  jet  is  cal- 


Tj  =  rnj  Vj 


where  mj  is  the  mass  flow  at  the  jet  is 

mj=pAjVj  * 


Using  Equation  (5)  and  the  integral  form  of  Equation  (10), 
the  jet  momentum  becomes 

X  =  1 


Tj -AttR J  [Pti-PjJ  xdx  • 

X=Xw  ^ 

Table  2  is  a  sample  of  the  data  reduction  sheets  which 
were  used  to  calculate  the  efficiency,  mass  flow,  and  jet 
momentum  for  this  program. 


2. 7. 2. 5  Augmentation  Ratio,  A 


The  augmentation  ratio  is  defined 
as  the  total  lift  divided  by  the  jet  momentum, 

A-  J=- 

A "  m  Vj  tt*) 

2. 7. 2. 6  Inlet  Ring  Thrust,  Ts 


The  inlet  ring  thrust  was  obtained 
by  a  numerical  integration  of  the  axial  component  of  the 
static  pressure  times  the  cross-sectional  area  of  the  inlet 
ring .  Thus , 

Rsfc 

Ts  =2tt J  Psrdr 

RS1  U3) 

where  Rs^  and  RS2  are  the  inner  and  outer  radii  of  the  inlet 

ring  respectively.  The  static  pressures  utilized  for  this 
calculation  were  the  average  of  pressures  measured  at  four 
azimuth  positions  on  the  ring. 

This  measurement  of  inlet  ring 
thrust  was  checked  by  determining  the  loads  on  the  three 
fixed-beam  assemblies  supporting  the  inlet  ring.  These 
data  were  found  to  agree  with  the  pressure  thrust  measure¬ 
ment  if  the  pressure  in  the  seal  between  the  duct  and  the 
inlet  ring  was  included.  Due  to  the  additional  complica¬ 
tion  of  evaluating  the  pressure  in  the  seal,  the  inlet  thrust 
measurements  presented  in  this  report  are  from  the  inlet 
pressure  measurements. 
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APPENDIX  II 


ANALYSIS  OF  THE  LOADS  AND  MOMENTS  ON  A  TWISTED  AND 
TAPERED  PROPELLER  BLADE  OPERATING  IN  A  DUCT  WITH 
AZIMUTHALLY  UNIFORM  FLOW 

This  appendix  presents  an  analysis  for  the  determin¬ 
ation  of  loads  and  moments  developed  in  a  twisted,  tapered, 
and  cambered  blade,  as  well  as  a  sample  calculation  for 
the  blade  used  in  this  program.  It  is  assumed  that  this 
analysis  is  performed  following  the  performance  analysis 
discussed  in  Section  1.2  so  that  the  distribution  of  the 
aerodynamic  loading  parameters  with  blade  radius  are 
available. 

The  forces  and  moments  on  the  blade  which  are  of 
principal  interest  are: 

1.  Centrifugal  Force 

The  three  orthogonal  components  of  the  force 
acting  on  the  propeller  blade  which  are  parallel  and  per¬ 
pendicular  to  the  axis  of  rotation  are  called  the  thrust, 
the  drag,  and  the  centrifugal  force.  The  thrust  load  is 
parallel  to  the  axis  of  rotation  and  may  be  calculated  by 
dividing  the  propeller  thrust,  Tp,  obtained  from  the  per¬ 
formance  calculations  by  the  number  of  blades.  This  thrust 
load  will  usually  be  small  for  GEM  propellers.  The  drag 
force  is  very  small  and  may  be  neglected  for  most  designs. 
Centrifugal  force  is  usually  by  far  the  largest  of  the 
forces  on  the  propeller  blade.  This  force  is  therefore  re¬ 
quired  to  determine  the  stress  on  the  blade  root  and 
attachment. 

2.  Out-of-Plane  Bending  Moment 


The  out-of-plane  moment  due  to  aerodynamic  loading 
can  be  balanced  by  coning  the  blades  in  the  direction  of  the 
thrust.  This  causes  a  moment  due  to  the  centrifugal  loading 


which  is  in  opposition  to  the  aerodynamic  moment.  Proper 
selection  of  the  coning  angle  will  result  in  a  significant 
reduction  of  blade  root  stresses. 

3.  In-Plane  Moments 

These  are  moments  on  the  blade  and  root  attach¬ 
ment  in  the  plane  of  rotation  due  to  the  torque  required 
and  possibly  due  to  a  component  of  the  centrifugal  force. 
Again,  these  moments  may  be  in  opposition  and  can  be  balanced 
to  reduce  blade  attachment  stresses. 

4.  Torsional  Moments 

Blade  torsional  moments  result  from: 

a.  Airfoil  section  pitching  moments. 

b.  Moment  due  to  the  location  of  the  airfoil 
section  aerodynamic  centers. 

c.  Moment  due  to  the  component  of  centrifugal 
force  which  results  from  blade  twist. 

(This  moment  is  sometimes  referred  to  as 
the  tennis  racket  or  the  dumbbell  effect.) 

By  proper  location  of  the  blade  section  aero¬ 
dynamic  centers,  the  torsional  moment  can  also  be  made  small. 

When  the  loads  and  moments  on  the  blade  and  the 
root  connection  are  established,  standard  stress  analysis 
procedures  may  be  applied  for  material  selection,  determin¬ 
ation  of  bolt  sizes,  retention  fitting  geometry,  etc.  These 
considerations  are  not  discussed  in  this  section. 

The  notation  used  in  this  section  is  illustrated 
in  Figures  26  and  27. 

1.0  ANALYSIS 

The  equations  required  to  determine  the  significant 
loads  and  moments  on  a  GEM  propeller  blade  are  presented 
in  the  following  subsections.  The  methods  used  for 
making  the  moments  small  are  also  given.  It  is  assumed 
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that  the  blades  are  rigid  for  this  analysis. 
1.1  Centrifugal  Force 

The  centrifugal  force,  Fop 

r 

x  =  1 

•(OR  }‘J  dm 


is  given  by 


(1) 


Fcf 


dr 


xdx 


(2) 


X=  Xta 


where  dm/dr  is  defined  as  the  change  in  blade  section  mass 
with  radius. 

1.2  Out-Of-Plane  Bending  Moments 

The  out-of-plane  bending  moment  is  due  to  the 
aerodynamic  loading  and  due  to  the  centrifugal  loading  on 
the  blades  which  are  inclined  at  the  coning  angle,  d0  . 
The  moment  due  to  aerodynamic  loading  is 

R 


M 


OA 


■/ 


r  d  T 


(3) 


f?H 


or,  in  a  form  more  suited  for  calculation, 

Moa=ra^M  r  j£LXdx 

b  J  dx 

The  moment  due  to  centrifugal  loading  is 

R 


(4) 


M 


oc 


dm  »  ~  2  — 


dr 


Q  rzdr 


(5) 


Rh 


44 


where  Z  is  the  distance  from  the  blade  section  centroids 
to  a  plane  which  is  perpendicular  to  the  axis  of  rotation 
and  through  the  reference  center  of  the  root  attachment. 

If  the  variation  of  Z  with  blade  radius  is  linear  and  the 
resulting  coning  angle  is  small, 


Z  =  rao 


Substituting  this  relation  into  Equation  (5)  and  writing 
in  terms  of  x, 


*  r 
Moc  =  (f2R)  RQo  J 


x  dx  • 


Equating  the  aerodynamic  moment  to  the  moment  due 
to  centrifugal  force,  for  zero  banding  at  the  root  attach¬ 
ment,  the  angle  a  can  be  found. 

c  p  f  7  xdx 

„  _  op  H  Jxh  d  x 


b  ['  H?X  dx 
Jxr a  dr 


If  the  variation  of  Z  with  r  cannot  be  approximated 
by  Equation  (6) ,  the  following  iterative  procedure  should 
be  followed: 


a.  Select  Z  at  the  tip  of  the  blade,  Zt 

b.  Locate  the  blade  section  centroid  at  each 

radial  station  for  the  required  pitch  setting 
and  twist  angle  in  reference  to  . 

c.  Determine  the  distance  from  the  blade  section 
centroid  to  the  plane  of  the  root  attachment. 

d.  The  moment  contribution  of  the  centrifugal 
force  on  each  blade  section  can  then  be 
determined.  If  the  total  flapping  moment 
due  to  centrifugal  force  is  found  to  differ 
from_the  aerodynamic  moment,  a  new  value 

of  Zt  assumed  and  the  calculation  pro¬ 

cess  is  repeated. 
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1.3  In-Plane  Moments 


The  in-plane  moments  on  the  blade  are  due  to 
propeller  torque  from  aerodynamic  loading  and  from  the  effect 
of  centrifugal  force  on  the  lateral  offset  of  the  blade 
section  centroids. 


1.3.1  In- Plane  Moment  Due  to  Torque 

The  moment  due  to  propeller  torque  acts 
opposite  to  the  direction  of  rotation  and  is  calculated 
as  follows: 


.  HP(55Q)  . 

Mia  -  ^ 


(9) 


1.3.2  In- Plane  Moment  Due  to  Centrifugal  Force 

The  component  of  the  centrifugal  force 
acting  in  the  plane  of  rotation  causes  a  moment  as  follows: 


X=XrA 


where  Yx  is  the  lateral  distance  from  the  blade  section 
centroids  to  the  vertical  projection  of  the  root  attach¬ 
ment  centroid. 

1.4  Airfoil  Pitching  Moment 

The  torsional  moments  on  the  blade  are  evaluated 
as  follows: 


1.4.1  Airfoil  Pitching  Moment 

The  pitching  moment  at  the  aerodynamic 
center  of  the  blade  sections  causes  a  blade  torsion  moment. 


M 


l '^/thrust"  b  J  tt*) 

where  y  is  the  lateral  distance  from  the  centroid  of  the 
root  attachment  to  the  25-percent  chord. 

1.4.2  Torsional  Moment  Due  to  Centrifugal  Force 
and  Blade  Twist 

Because  of  twist  and  taper,  the  centrifugal 
force  acting  through  the  centroids  of  each  blade  section 
generates  a  torsional  moment. 

The  centrifugal  force  of  element  (dydr)  is 


dRF  =  dfidxdrhyft2(r**y*)  * 


where  de  •  Mass  Density  of  Material. 


ial  direction  is 


The  component  of  that  force  in  the  tangent- 


dFcr  =  ^dedydrhyftz(rz+yz) 2  j  sin(J)) 


The  torsional  moment  is  the  double  integral 
of  the  sine  component  of  this  tangential  force  across  the 
limits  of  the  blade  leading  to  trailing  edges  and  root  to 
tip  respectively: 

R  Vie 

Mrc  =  /  /  ysin^sina  [dedydrhyQ^r^y  )*j-  ^ 

rvA  Yte 


48 


Substituting  for 
collecting  terms, 


sin0y 


y 

(r*+y*)* 


into  Equation  (17)  and 


X  =1  Yle 


Mrc  =deQ.*  rJ~  J hyysin0*dydx 


x=XrA  Yte 


(18) 


After  solving  for  the  loads  and  moments  as  indicated,  one  may 
apply  standard  stress  analysis  procedures  for  material 
selection,  bolt  sizes,  retention  fitting  geometry,  etc. 


2.0  SAMPLE  CALCULATIONS 


This  section  presents  the  stress  calculations  used  in 
the  design  of  the  optimum  propeller. 


The  design  conditions  of  the  optimum  blade  are  presented 
in  Table  3.  The  blade  chord  and  twist  are  given  in  Figure 
25. 

The  blade  tip  speed, Q  R,  varies  with  nozzle  height  and 
was  obtained  from  performance  calculations  for  18  p.s.f. 
disc  loading  and 


QR  = 


Ct  p 


(19) 


For  this  analysis,  however,  we  will  consider  a  design  con¬ 
dition  of  h/d  =  0.27.  This  is  the  maximum  height  at  which 
the  test  unit  will  operate  for  any  length  of  time.  The 
design  values  selected  were  for  an  h/d  of  0.27  and 
Qr  -  552  f.p.s.,  HP  =  14.9  horsepower,  D.L.  =  18  p.s.f. 

It  may  be  seen  from  Figure  28  that  the  blade  can  be 
conveniently  analyzed  in  two  segments: 
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Segment  I  -  Twist  and  Taper  Airfoil  Section  from 
x=xa  to  x=i . 

Segment  IX  -  Tapered  Root  Section  from  X  =  Xra  to  X  =  XH 

Numerical  integrations  were  used  in  the  solution  of  equations 
For  this  purpose,  Segment  I  was  divided  into  five  equal  sect¬ 
ions  as  shown  in  Figure  28. 

Figure  28  also  indicates  the  sections  of  Segment  II. 

This  segment  is  treated  as  a  rectangular  prism  and  two 
triangular  pyramids.  The  centroids  of  these  sections  are 
located  at  radial  locations  Xa  and  Xb  respectively. 


2.1  Centrifugal  Force 


The  numerical  integration  of  Equation  (1)  is 

X=1 

Am  \ 


CF 


A  r 


A  x 


A  m 


X=XrA 


(20) 


By  definition,  =  (Mass  Density)  (Cross-Sectional 

Area),  slugs/foot. 


2.1.1  Centrifugal  Force  for  Segment  I 


The  cross-sectional  area,  S*,  for  a  given 
airfoil  is  proportional  to  the  square  of  the  chord.  There¬ 
fore,  for  a  blade  with  constant  airfoil  section,  the  area 
at  any  station  X  may  be  written  as 


Sx  =  ks  Cx 


The  area  of  a  16-709  airfoil  section  with  a  10-inch  chord 
was  determined  with  a  planimeter  to  be  7.0  square  inches. 
Therefore, 

S* 


ks  = 


Cx2 

=  0.07 
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With  the  chord  values  of  Figure  25,  the  area  of  airfoil 
sections  at  each  station  was  calculated,  and  the  results 
are  summarized  in  Table  4. 


From  Reference  6,  the  density  of  white  ash 
42  pounds/cubic  foot 


Am 


42 


x  32.2 
Am 


A  r 

The  calculated  values  for 
presented  in  Table  4. 


A  r 


Sx|=  1.304 Ax 

at  each  station  are  also 


1 


From  Table  4,  ^  ^r~r  XAX 


X=Xw 


Ar 


0.00264. 


2.1.2  Centrifugal  Force  for  Segment  II 

The  following  values  are  needed  for  the 
rectangular  prism  portion: 

X=  1734°  °'282 
Ax=  -  =0.352 


17.34 


Am 

Ar 


=1.304  S*  =1.304 


(0.88 

(3.0) 

12 

12  1 

=  0.0119 


slug 

ft. 


For  the  triangular  pyramid  portions  of  Segment  II, 

x  =  iff*  =0-39 

AX=lf§4°a19 


Ar 


2 

0.88\ 

2  1.031 

3 

12 

3  12  ) 

ft- 
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**?»<**•  *-  «•  "*■ 


for  Segment 


The  summation  of  XAX 


Am 
A  r 


II  is  X=Xh 

2_xAx  M  =0.282(0.352)0.0119  ♦2(0.39)(0.19)(0.0035) 

X=Xr* 

=  0.00116  ♦  0.00052 
=  0.00168  • 

The  total  centrifugal  force  may  now  be  calculated  by 
substituting  the  above  values  into  Equation  (20) : 

Fcf  =(552)'  jo.00264  *0.00168] 

=  1320  lb. 

2.2  Out-of-Plane  Moment  Due  to  Aerodynamic  Loading 
Equation  (4)  may  be  written  as 

pM'  X=1 


Mo» 

obtaining 


_  RSpI 

fc 

ACt 

Ax 


V  ACt 
L—  Ax 


xAx 


X=XH 


from  performance  calculations. 
x-1 


From  Table  4,  ^  XAXr  0.0187 


X«*. 


Mo*  =  1445,(5.21)(0.00238)(552)g(00187) 


=  25.5  ft.- lb.  • 
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2.3  Blade  Coning  Angle  for  Zero  Out-of-Plane  Bendin* 


Assuming  the  coning  angle  is  small  and  rewriting 
Equation  (8) , 

X=1 


-^pI^Txax 

U  X«Xm 


?AX<X‘>  If  -£x'Ax%f 


From  Table  4, 


t  =  0.0187 

Ax 

X  =  1 

^x'Ax^p-  =0.001676  • 

X«Xh 

From  centrifugal  force  calculations, 

K=Xh 

Y_  A  xx*  =0.0119  •  0. 3 5(0.28)' *  2(0.0035X0.1 9) 0. 39 ’ 


C=Xr» 


=  0.00035  *  0.000203 
=  0.000528 

■?(5.21)0.00238(0.0187) 
0.001676  ♦  0.000528 

=  0.263  rad. 


=  1.51  degrees 
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With  the  above  value  for  Oo  ,  the  distance,  Zx  , 
was  obtained  by  using  the  method  described  in  Section  1.2  of 
this  appendix.  The  resulting  moment  due  to  centrifugal 
force  was  found  to  be  in  close  agreement  with  the  aero¬ 
dynamic  moment. 

2.4  In- Plane  Moment  Due  to  Torque 

As  mentioned  previously,  the  design  maximum  power 
required  for  this  propeller,  obtained  from  performance  cal¬ 
culations,  is  14.9  horsepower  at  a  tip  speed  of  552  f.p.s. 

Substituting  these  values  into  Equation  (9) , 

k,  HP(550) 

M"— 5 - 

14.9(550) 

SS2 
144  S 

=  20.9ft:lb.  or  per  blade  =5.24ftrlb. 

Since  the  rotation  is  counterclockwise,  the  moment  due  to 
torque  acts  clockwise. 

2.5  In-Plane  Moment  Due  to  Centrifugal  Force 

The  only  contribution  to  the  in-plane  moment  is 
from  Segment  I,  since  Segment  II  is  symmetrical  about  the 
root  attachment  centroid.  The  distance*  yx  ,  was 

obtained  by  first  determining  the  centroid  of  each  blade 
section  by  the  standard  area  moment  method.  The  centroids 
for  each  blade  section  were  found  to  lie  on  the  airfoil 
mean  line  at  the  50-percent  chord.  Values  for  y*  were 

then  measured  from  a  layout  of  the  blade  and  are  presented 
in  Table  4. 


The  contribution  of  each  blade  section  can  now  be 
determined.  Combining  Equations  (2)  and  (10)  in  a  form  for 
numerical  integration, 


5A 


r  v*.  yj  -  “ i1  '■•I'jf . 


From  Table  4, 

X=  1 

V  x  Ax^x  =  °-0000245 

X  =  XH 

for  Q  R  «  552  ft. /sec. 

Mk  =(552)'[0.0000245] 

=  7.45  f t-tb. 

This  moment  acts  counterclockwise,  since  centroids  lie 
ahead  of  root  attachment  centroid  projection  (see  Figure  26). 


The  residual  in-plane  moment  is 


Mlq 

=  5.24 

ft- lb 

cw 

Mic 

=  7.45 

ft- lb 

ccw 

Mr 

=  2.21 

ft-lb 

ccw 

2.6  Torsional  Moment 


2.6.1  Torsional  Moments  Due  to  Airfoil  Pitching 
Moment 


Airfoil  section  pitching  moment  coefficients 
were  obtained  from  Reference  6.  The  ratio  of  induced  velocity 
to  tip  speed  and  the  local  solidity  Ox  were  obtained  from 
performance  calculations.  These  data  are  included  in  Table  4. 
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Equation  (13)  may  be  written  as 


p/cm,j^)%xs]Ax 


X=Xh 


From  Table  4,  Ox  Cm,  +X  AX=  -0.00180 

X=X« 

MtC  _(552)* 0.00238 7T qqq^qI 
4  2(4)  L  J 

=  “1.22  ftrlb.  (Nose  Down). 

2.6.2  Torsional  Moment  Due  to  Aerodynamic  Lift 

For  this  blade,  the  centroid  of  the  root 
attachment  is  not  on  the  25-percent  chord;  hence,  a 
torsional  moment  arises  from  the  aerodynamic  lift. 

From  blade  layout,  y  =*1.27  inches 
(scaled  dimension) . 

From  Equation  (14) , 

(Mt)thrust  =  qj~ '  y 

_  93.8  v  1.27 

"  4  *  12 

=  2.49  ftrlb.  (Nose  Up). 


2.6.3  Torsional  Moments  Due  to  Centrifugal  Force 


2. 6. 3.1  Torsional  Moment  Due  to  Segment  I 
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Since  both  airfoil  thickness  and 
length  are  functions  of  chord,  the  values  of  hy  and  y  from 
Figure  27  can  be  written  as  percent  of  chord,  and  hence 
hy  and  y  remain  constant  at  any  x  value. 

For  numerical  integration  purposes, 
the  airfoil  section  is  assumed  to  be  divided  into  five 
segments,  each  equal  to  20-percent  of  the  chord. 


The  root  attachment  centroid 
projection  lies  along  the  50-percent-of-chord  line  at 
x  *  xjj.  values  of  y  and  hy  are  obtained  from  the  airfoil 
ordinates  and  are  also  presented  in  the  table  below. 


STA 

y 

c 

hi 

C 

yh 

C3 

1 

0.4 

0.03186 

0.0051 

2 

0.20 

0.08126 

0.0032 

3 

0 

0.0900 

0 

4 

-0.20 

0.07904 

0.0031 

5 

-0.40 

0.03776 

0.0060 

-4^  =0.20 

y=&  c 

Y_  y*hyAy  =  0.00348  CS 

V  =  +e 

The  numerical  integration  form  of 
Equation  (18)  for  Segment  I  is 

y  =ie  X= 1 

=d*azRVy*hyAyy"  sin  ©x  Ax 
y  =  te  X=Xh 

x=1 

0.00348^~C*4  sin©*Axj  . 

X=Xh 


Mtc 


Mtc  =  deftzR 
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-r-  r  ^ 


X  =  1 

From  Table  4,  Yc.  sina  Ax  = -0.000986 
X  =  XH 

Mrc  =1,304(382)*1.445(0.00348)(-0.000986) 

=  -0.94ftrlb.  (Nose  Down). 

2 . 6 . 3 . 2  Torsional  Moment  Due  to  Segment  II 

The  numerical  equation  form  of 
Equation  (18)  for  Segment  II  is 

Mrc=de^zR/_  /_  hyyAyAxsin0x 

y*te  x«x« 

Q^for  this  segment  is  constant 
Sin  0x  =  -0.301. 

Using  the  values  for  the  rectangular  prism  and  the  triangular 
pyramid  sections  described  in  Section  2.1.  and  shown  in 
Figure  32,  the  torsional  moment  is 

I  n  In  7R12m  e.  \l 

0.46-0.11 


Mrc  =  1. 304(382 H1.445)(-0.301)  2  \ ^ 


co 

CO 

o 

[0.75f 

1.5 

1  12 

12  1 

12 

*1 


0.88)11.83 |2  2(1.012 


12 


12"  3  if)  3 10'46 -°'27! 


=  -8280  [0.000025  +  0.00  0016 
=-0.34ftrlb.  (Nose  DowtO- 

2.6.4  Total  Torsional  Moments 


Summing  up  torsional  moments  due  to  centri¬ 
fugal  force  of  Segments  I  and  II, 
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mTC  *  1.28  foot-pounds  (nose  down). 


The  total  torsional  moment  acting  on  the 
root  attachment  fitting  is  the  summation  of 


Mtc  =Mrci.  *MrTumsT  *  Mtc 

=-1.22*2.49  -1.28 

Mr  s  0  ft-  lb.  • 


3.0  BLADE  DESIGN  LOADS 


In  summation,  the  following 

Centrifugal  Force 
Out-of-Plane  Moments 
Coning  Angle, d0 
In- Plane  Moment 
Torsional  Moment 


results  were  obtained: 

1320  pounds 

0  foot-pounds 
1  degree  30.5  minutes 
13.5  foot-pounds 
0  foot-pounds 


These  loads  were  used  to  design  the  root  attachment  and 
hub  of  the  optimized  propeller  tested  in  this  program. 
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APPENDIX  III 


PROPELLER  OPTIMIZATION  PROCESS 
AS  APPLIED  IN  THIS  PROGRAM  TO  THE  OPTIMIZED  PROPELLER 


The  following  development  has  a  twofold  purpose: 
first,  the  design  considerations  and  procedures  for  the 
propeller  used  in  this  program  will  be  given;  second, 
the  procedure  of  Section  2.1  will  be  followed.  Thus,  this 
development  can  be  used  as  an  illustrative  example. 

It  is  noted  that  in  some  cases  the  design  procedure 
used  deviated  from  the  design  procedure  of  Section  2.1. 

This  occurred  since  the  design  of  the  optimized  propeller 
followed  the  work  of  Reference  1. 

The  diameter  of  the  optimized  propeller  was  made  to  fit 
the  dimensions  of  the  test  duct  facility.  The  local  solidity 
of  the  optimum  propeller  was  made  the  same  as  the  straight- 
bladed  propeller  at  the  0.8-blade  radius  for  the  purpose  of 
comparison . 

As  a  result  of  the  work  of  Reference  1,  there  was  an 
accumulation  of  propeller  performance  and  GEM  performance 
data  obtained  at  various  test  conditions.  In  particular, 
some  attention  was  focused  at  the  case  where  h/te  =  i.oo 
which  represented  an  h/d  =  0.27,  and  —  .TL.d.b —  =  1.37 

for  the  case  of  Sj/s  =  0.80.  Thus,  the  newly  designed 
propeller  was  optimized  for  operation  at  this  height.  With 
a  reduced  pitch  setting,  this  propeller  also  operates  at 
near  optimum  performance  at  lower  heights. 

Step  1:  Optimum  Airfoil  Section  Performance 

It  is  necessary  to  select  an  airfoil  with  the 
best  aerodynamic  performance  and  to  prepare  the  character¬ 
istics  of  the  selected  airfoil  section  in  a  form  which  can 
be  used  in  the  calculations.  The  procedures  used  to  achieve 
these  results  are  discussed  below. 
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The  aerodynamic  characteristics  of  various 

airfoils  were  summarized,  and  a  plot  of  -~LJ_  versus 

r  Cd 

section  lift  coefficient  for  these  airfoils  is  shown  in 

Figure  29. 

c  — 

The  maximum  value  of  — is  attained 

C  c  - 

by  the  airfoil  16-709,  and  the  values  of  — ■  on  each 

Cs  ~  w 

side  of  the  maximum  value  ~---a  -  are  seen  to  be  higher 

than  that  of  the  other  airfoils  investigated.  Thus,  this 
airfoil  would  give  the  best  efficiency  of  the  airfoils 
investigated,  as  explained  in  Section  2. 


The  lower  surface  of  the  16-709  airfoil  is 
almost  a  straight  line,  which  makes  its  fabrication  simpler 
and  less  costly.  This  airfoil  causes  an  aerodynamic 
pitching  moment  which  must  be  accounted  for  in  the  stress 
calculations  as  explained  in  Appendix  II. 


The  airfoil  characteristics  for  the  16-709 
airfoil  are  presented  in  Figure  30.  These  data  are  presented 

in  the  form  and  plotted  against  angle  of  attach 

for  various  Mach  numbers  in  Figure  31.  The  angle  of  attack 

at  which  -  is  a  maximum  was  read  from  the  curves 

Cd 

of  Figure  31  and  was  plotted  against  Mach  number  in  Figure 
32.  This  curve  is  used  in  the  calculations. 


Step  2:  Propeller  Tip  Sneed 

As  discussed  in  Section  1.2.2  of  this  report, 
a  desired  rotational  speed  is  generally  available.  The 
nondimensional  propeller  performance  parameters  generally 
do  not  vary  significantly  for  propeller  tip  speeds  from 
400  to  650  f.p.s.,  and  therefore  calculations  based  on  any 
speed  in  this  range  will  be  almost  identical.  A  propeller 
tip  speed  of  460  f.p.s.  was  used  for  the  sample  calculation. 
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The  propeller  for  this  program  was  optimized 
for  a  nozzle- area- to-propeller- area  ratio,  Sj/Sp,  of  unity. 
The  propeller  was  assumed  to  be  divided  into  the  radial 
increments  shown  in  Figure  28.  The  radial  stations  at  the 
middle  of  these  increments  are  listed  in  Table  5. 


Step  4:  Axial  Velocity  Ratio 

As  was  previously  mentioned,  the  propeller 
solidity  of  the  optimized  propeller  is  equal  to  the  non- 
twisted  propeller  solidity  at  the  0.8  blade  radius.  Also, 
the  optimized  propeller  is  to  produce  uniform  axial 
velocity,  and  therefore  the  axial  velocity  ratio  will 
be  constant  with  radius  for  this  propeller.  Thus, 

Equation  (45)  of  Section  2.2  can  be  solved  for  the  axial 
velocity  ratio  for  the  optimized  propeller  if  the  follow¬ 
ing  assumptions  are  made: 


v 


a .  The  value  of  1  /  1  + 

independent  of  Vp/ Q  R.  That 


is  practically 


is, 


V1  =  1 


b.  The  product  of  the  drag  coefficient  and 
the  axial  velocity  ratio  is  very  small 
as  compared  to  the  section  lift  coeffic¬ 
ient.  That  is, 

Cd(oRx)  ^  c< ' 

c.  Inlet  loss  factor  is  zero.  (No  data  on 
this  factor  were  available  at  the  time 
of  this  analysis.  See  Figure  9  for  a 
better  value  for  this  factor.)  With 
these  assumptions, 


Vp  . 

or  ‘ 


UX  v  IS 

t  Qir 


To  solve  this  equation,  a  value  for  Ts/Tp 
must  be  established.  It  was  assumed  that  there  would  not 
be  a  significant  influence  of  the  propeller  design  on  this 
parameter,  and  therefore  the  data  obtained  with  the 
straight-bladed  propeller  shown  in  Figure  9  could  be  used 
for  the  optimized  propeller  design.  As  shown  in  this 
figure,  this  assumption  was  substantiated  by  the  later 
testing.  A  value  for  Ts/Tp  of  0.23  was  used  for  the 
design . 


The  value  for  C$  to  be  used  in  the  equation 
for  Vp/  Q  R  is  the  Cj|  obtained  at  the  optimum  angle  of 
attack.  It  is  first  necessary  to  calculate  the  section 
Mach  number  as 


m  -  XQR 
Mx~  1120 


_  0.8(460) 
1120 

=  0.32 


Then  from  Figure  32  the  optimum  angle  of  attack  of  the 
16-709  airfoil  is  3.2  degrees.  From  Figure  30  the  lift 
coefficient  at  this  angle  of  attack  is  seen  to  be  0.81. 


as 


The  axial  velocity  ratio  can  now  be  calculated 


Vp  _ 
OR  " 


Ox  v  Ts  \  z 

ITxQ  17} 

(^(0.8X0.81)0.23 


L 

z 


=  0.127  . 
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This  axial  velocity  is  kept  constant  at  all  radial  stations 
as  shown  in  column  3  of  Table  3. 


Step  5:  Blade  Section  Mach  Number 


The  Mach  number  at  each  radial  station  can  be 
calculated  from  the  given  tip  speed  of  460  f.p.s.  as 

*^60) 

1120 

and  placed  in  column  4  of  Table  5. 

Step  6:  Optimum  Angle  of  Attack 

The  optimum  angle  of  attack  for  each  station 
can  be  read  from  Figure  32  for  the  calculated  blade  station 
Mach  numbers  and  placed  in  column  6  of  Table  5. 


Step  7;  Propeller  Inflow  Angle 
Propeller  inflow  angle, 


can  be  calculated 


as 


(j)  =  tan" 


Vp 


XHR  J 

and  placed  in  column  7  of  Table  5. 

Step  8;  Propeller  Blade  Pitch 

The  blade  pitch  can  be  obtained  from 

0*  =  CXopr  +  0  • 


These  data  are  shown  in  column  8  of  Table  5. 

Step  9:  Airfoil  Lift  and  Drag  Coefficients 

and  Cd  can  be  obtained  from  Figure  30 
for  the  angles  of  attack  of  colunn  6  of  Table  5  at  the 
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Mach  numbers  of  column  4.  Columns  9  and  10  of  this  table 
present  Cj|  and  Cd  for  these  conditions. 


Step  10:  Solidity 

For  the  case  of  Sj  »  Sp  and  k  »  0,  Equation 
(45)  reduces  to  2 

VP 


Gx_  IflR  1 

The  values  obtained  from  previous  steps  are  substituted  into 
this  equation  and  the  results  are  shown  in  column  15  of 
Table  5. 

Step  11:  Propeller  Thrust  and  Power  Coefficient 

The  columns  of  Table  5  are  set  up  according 
to  Equation  (23) .  Column  20  gives  the  incremental  thrust 
coefficient.  Adding  this  column  and  multiplying  by  the 
incremental  span, 


T”Z_  Ax 


Ct  -  (0.324)  (0.108)  -  0.0350. 


The  thrust  produced  for  this  thrust  coefficient  is 


Tp  =  CrPSP(nR)Z 
=0.035(0.00238X5.21)460* 
=  91.6  lb.. 


65 


Propeller  power  is  calculated  similarly  to  be 


AX 


=0.0454(0.108) 


=0.0049 . 


The  horsepower  required  for  this  coefficient  is 


SHP  = 


CpPM  Sp 

550 

0.0049(0.00238)4603(5.21) 

550 


=  107  hp  • 

As  a  result  of  this  analysis,  the  propeller  characteristics 
are  designated.  The  characteristics  of  the  optimized 
propeller  are  shown  in  Figure  25. 
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TABLE  1 


GEM  DUCT- TEST  GEOMETRIC  PARAMETERS 


Total  Base  plus  Jet 

Cross-Sectional  Area,  S 

6.6 

Square 

feet 

Inlet  Cross-Sectional 

Flow  Area,  Si 

5.21 

Square 

feet 

Propeller  Plane 

Cross-Sectional  Area,  Sp 

5.21 

Square 

feet 

Exit  Nozzle  Area,  Sj 

Configuration  1 

5.21 

Square 

feet 

Configuration  2 

3.90 

Square 

feet 

Configuration  3 

2.00 

Square 

feet 

Configuration  4 

0.65 

Square 

feet 

Propeller  Hub  Radius  to 

Propeller  Tip  Radius  Ratio, 

Rh 

"r 

0.46 

Length  from  Propeller  Plane 
to  Duct  Exit  Plane,  Z 

4.42 

feet 

TABLE  2 


TYPICAL  DATA  REDUCTION  SHEET 


DATE  20  JUL  62 
TEST  NO.  97 
D.L.  »  16 


h/d  =0.15  rj j  =62.7  7 o 

2/p  =900  Tj  =  26.50  ib 
IHP  =8.70  F  =  238.5  crs. 


Ex 

\ 

Ef 

* 

©»© 

(D  X  © 

mm 

.0198. 

.  507 

o.  10 

e.  120 

0.16,2 

2.94-8 

o.eio 

3.  04-0 

0.24-8 

4-.  224- 

0,210 

V 

CD 

O 

ptj‘PJ  Ptj 


ax  ®»®  4®  V© 


©  ®  ®  ©  ®  ®  ®  ®  ® 


t.io  (1.50 


0.04  0.018  ^°  3U5  o .zo8 


2.10 


3.oo  1 2.T 


3.  )0  |2. 


3.0  12.3 


0.50 


0.60 


0.70 


0.80 


0.90 


Looid  i&.P*) 


.024-1®  Z-UZltfL®  13.12KB) 


INTERNAL  EFFICIENCY,  V\  =  Y.  ®] 

VOLUME  FLOW,  F  =  27^^© 

MOMENTUM,  Tj  =  47rRzY® 


TABLE  3 


OPTIMIZED  PROPELLER  CHARACTERISTICS 


Number  of  blades,  b 

4 

Outer  radius,  R 

1.445  feet 

Hub  Radius,  R^ 

0 .665  feet 

Centroid  root  attachment  fitting,  Ra 

0.159  feet 

Blade  chord,  c 

function  of  radius 

Blade  pitch,  0b,  (at  X  ^  0.46) 

17.50°,  19.4° 

Blade  twist,  @b, 

function  of  radius 

Maximum  design  propeller  thrust,  Tp 

93.8  pounds 

Maximum  propeller  disc  loading  Tp/Sp 

18  p.s .f . 

Design  tip  speed,  ft  R 

460  f.p.s. 

Planform 

tapered 

Airfoil 

NACA  16-709  series 

BLADE  PARAMETERS  FOR  DETERMINATION  OF  BLADE  STRESSES 


-0.00519  -0.00215  -0.00101  -0.00050  -0.00028  -0.000986 


DIGITAL  PROCESS  FOR  THE  DETERMINATION  OF  PROPELLER  THRUST  AND  POWER  COEFFICIENT 


to 

r-4 

00 

CO 

CM 

r-4 

rH 

O 

O 

O 

o 

o 

o 

• 

• 

• 

• 

• 

rH 

t-4 

rH 

rH 

rH 

to 

CM 

at 

o 

© 

CM 

00 

CO 

rH 

rH 

rH 

O 

CM 

00 

to 

m 

tO 

CO 

CM 

rH 

o 

<r 

to 

to 

00 

o 

o 

O 

O 

o 

-a 

v  v 

o 

o 

o 

O 

o 

• 

o 

• 

o 

• 

o 

• 

o 

• 

o 

© 

w  * 
© 

• 

o 

• 

o 

• 

o 

• 

o 

• 

o 

CM 

to 

o 

to 

00 

oo 

o 

o 

to 

o 

o 

o\ 

at 

00 

/ 

to 

at 

r-4 

CO 

rH 

rH 

o 

o 

o 

V. 

V 

tO 

to 

r- 

at 

o 

o 

o 

o 

o 

o 

< 

o 

o 

O 

o 

rH 

• 

o 

• 

o 

• 

o 

• 

o 

• 

o 

©  © 

• 

o 

• 

o 

• 

o 

• 

o 

• 

o 

o 

o 

o 

o 

oo 

o 

o 

o 

00 

rH 

rH 

rH 

rH 

o 

r-4 

r-4 

rH 

o 

00 

00 

00 

oo 

00 

^ao\ 

00 

00 

00 

oo 

• 

• 

• 

• 

• 

• 

• 

• 

• 

o 

o 

o 

o 

o 

Vc!7  © 

o 

o 

o 

o 

i 

Ot 

00 

00 

CM 

tO 

O 

CM 

00 

tO 

• 

• 

• 

• 

* 

tO 

o 

r^. 

m 

CO 

rH 

o 

CM 

CM 

rH 

rH 

rH 

rH 

rH 

O 

O 

* 

o 

o 

Ot 

©  © 

• 

o 

• 

o 

1 

1 

i 

CO 

CO 

to 

00 

CM 

• 

• 

00 

tO 

M3 

CO 

rH 

• 

• 

• 

r-4 

rH 

CTt 

oo 

O'. 

00 

00 

tO 

to 

A/ 

to 

r-l 

o- 

CM 

N 

CL 

CO 

at 

i". 

at 

co 

to 

O 

o 

o 

X 

CO 

-a 

tO 

00 

rH 

to 

•sT 

CM 

o 

i 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

a 

o 

o 

o 

rH 

CO 

CO 

CO 

CO 

CO 

tO 

o 

CO 

to 

o 

® 

CO 

to 

to 

CM 

rH 

rH 

<}• 

CM 

tO 

CO 

rH 

O 

at 

<1- 

o 

r» 

to 

CO 

rH 

rH 

rH 

rH 

o 

CM 

• 

CM 

• 

rH 

• 

rH 

• 

rH 

• 

@t5 

® 

• 

o 

• 

o 

• 

o 

• 

o 

• 

o 

O 

O 

o 

o 

o 

to 

at 

oo 

CM 

to 

o 

to 

at 

£ 

£ 

at 

rH 

CM 

CM 

CO 

CO 

CO 

1 

o 

rH 

• 

O 

• 

O 

• 

o 

• 

o 

• 

o 

©<§ 

• 

o 

• 

o 

i 

i 

i 

f>- 

f'. 

I-'- 

n- 

( 

5) 

at 

CM 

CM 

CM 

CM 

CM 

«< 

CM 

rH 

rH 

• 

o 

rH 

• 

o 

rH 

• 

o 

rH 

• 

o 

rH 

• 

o 

(9)  ® 
^  © 

• 

o 

to 

• 

o 

i 

i 

i 

Ml- 

CM 

o 

00 

to 

rH 

CM 

CO 

CO 

<fr 

m 

tO 

00 

at 

2 

1-4 

w 

H 

> 

> 

• 

• 

• 

• 

• 

© 

H 

H 

H 

o 

o 

o 

o 

o 

M 

</) 

H 

t— 1 

M 

> 

> 

H 

H 

H 

H 

71 


ISOLATED  PROPELLER  EFFICIENCY  ,  '/p,  PERCENT 


FIGURE  4:  PROPELLER  THRUST  DATA  FOR  VARIOUS  Sj/S  CONFIGURATIONS ,  OPTIMIZED  PROPELLER 


SOO  0 


FIGURE  5:  PROPELLER  POWER  DATA  FOR  VARIOUS  Sj/S  CONFIGURATIONS,  OPTIMIZED  PROPELLER 


STRAIGHT  BLADED 


GEM  CLEARANCE  ~~  NOZZLE  AREA  RATIO,  TTdh/Sj 


FIGURE  7:  DUCT  LOSS  FACTOR  FOR  VARIOUS  JET  AREA  CONFIGURATIONS 


P/  =  I6PSF. 
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THEORY  OF 
I  REFERENCE 


FIGURE  15:  COMPARISON  OF  PREDICTED  PERFORMANCE  WITH  PERFORMANCE 


INLET  LOSS  THEORY 
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FIGURE  19:  COMPARISON  OF  PREDICTED  AND  ESTIMATED  ANGLE  OF  ATTACK  OF  ROOT 

SECTION  (X  of  0.51),  OPTIMIZED  PROPELLER 


h/d  s  0.15  ' 

INLET  LOSS  THEOR 
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H  4. 5 -INCH  INLET  RING  AND  OPTIMIZED  PROPELLER 


LOCAL  BLADE  PITCH  ANGLE  ,  Ox,  DEGREES 


f*, 

I  FIGURE  24:  BLADE  CHARACTERISTICS  OF  STRAIGHT  BLADED  PROPELLER 
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LOCAL  BLADE  CHORD,  Cx  ,  INCHES 


LOCAL  BLAOE  PITCH  ANGLE,  Ox,  DEGREES 


REF.  PLANE  CONTAINING 
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PROPELLER  NOTATION 
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FIGURE  27:  NOTATION  FOR  THE  DETERMINATION  OF  PROPELLER 

TORSIONAL  MOMENTS 
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AIRFOIL  SECTION  EFFICIENCY  PARAMETER  Cx/C<j 


AIRFOILS  SHOWN  ARE  NACA  SECTIONS 

EXCEPT  WHERE  NOTED. 
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SECTION  LIFT  COEFFICIENT,  C^ 


FIGURE  29:  AIRFOIL  SECTION  EFFICIENCY  OF  SEVERAL  AIRFOILS 


M=  0.2  ( •  st.) 


U  0.2  0*  0.6  0.8  1.0 


SECTION  MACH  NUMBER,  M 

FIGURE  32:  VARIATION  OF  THE  OPTIMUM  ANGLE  OF  ATTACK  OF  THE 

16709  SERIES  AIRFOIL  WITH  MACH  NUMBER 
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